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Abstract 

We investigate the signatures of neutrinos produced in the annihilation of WIMP dark matter 
in the Earth, the Sun and at the Galactic centre within the framework of the Inert Doublet 
Model and extensions. We consider a dark matter candidate, that we take to be one of the 
neutral components of an extra Higgs doublet, in three distinct mass ranges, which have all been 
shown previously to be consistent with both WMAP abundance and direct detection experiments 
exclusion limits. Specifically, we consider a light WIMP with mass between 4 and 8 GeV (low), 
a WIMP with mass around 60-70 GeV (middle) and a heavy WIMP with mass above 500 GeV 
(high). In the first case, we show that capture in the Sun may be constrained using Super- 
Kamiokande data. In the last two cases, we argue that indirect detection through neutrinos is 
challenging but not altogether excluded. For middle masses, we try to make the most benefit 
of the proximity of the so-called 'iron resonance' that might enhance the capture of the dark 
matter candidate by the Earth. The signal from the Earth is further enhanced if light right- 
handed Major ana neutrinos are introduced, in which case the scalar dark matter candidate 
may annihilate into pairs of mono-energetic neutrinos. In the case of high masses, detection of 
neutrinos from the Galactic centre might be possible, provided the dark matter abundance is 
substantially boosted. 

1 Introduction 

Cosmological observations concur to indicate that about 25 % of the energy density of the universe 
is made of dark matter (DM) [HE]. The simplest and most popular paradigm postulates that dark 
matter is a weakly interacting massive particle or WIMP [3^ 4J. The reason is twofold. Firstly, a 
stable particle with annihilation cross section in the picobarn range, characteristic of weak interac- 
tions, would have a relic density in agreement with observations, ^dm ^ ^/{^^) ^ 0.25. Secondly, 
supersymmetric extensions of the Standard Model predict the existence of a dark matter candidate 
which, generically, is a spin 1/2 particle with weak scale interactions, the neutralino. There is also 
a popular spin 1 WIMP candidate, in the form of the lightest Kaluza-Klein partner of the photon 
in models with Universal Extra Dimensions (UED) [5j. To complete the list, one may consider the 
Inert Doublet Model (IDM), a very simple extension of the Standard Model introduced in ref. [6]. 
This model has two Higgs doublets and, to prevent FCNC, a discrete Z2 symmetry. The dark 



matter candidate is the lighest neutral component of the extra Higgs doublet Clearly, it is not the 
ambition of the IDM to compete with, say, the MSSM or UED models. This model is however very 
simple and, nevertheless, it does have interesting phenomenological consequences, some of which 
have been emphasized in recent works, refs. [SIEIITO]. Furthermore, its dark matter features are not 
less interesting, both from the perspective of direct detection and indirect detection, be it through 
gamma rays fTTj I12|. neutrinos [13j, or antimatter [14j. 

To complete the analysis of IDM dark matter phenomenology initiated in refs. [HI EJ [TOl [HI [121 
[T3t HH [m HE], we discuss further the signatures of the IDM dark matter candidate in neutrinos]^ 
Specifically, we consider three distinct dark matter mass ranges which we call in the sequel, the 
low, middle and high mass ranges. All are consistent with WMAP abundance and we consider the 
potential for indirect detection of these IDM dark matter candidates through the flux of neutrinos 
that might be produced in its annihilation in the Sun (low mass range), in the Earth (middle mass 
range), or at the Galactic centre (GC) (high mass range). 

As discussed in ref. [16j, the IDM dark matter candidate may be rather light, ttidm ^ few GeV 
(low mass range). Because of its potentially large cross section on nuclei, it has been advocated 
as one of the possible explanations for the recent combined DAMA/Nal and DAMA/Libra results, 
ref. [16]. Such a candidate may be a source of gammas from the GC fT6( [T8] and neutrinos from 
the Sun [191 [20]. Here, we discuss the constraint posed on the IDM candidate by the limits from 
the Super-Kamiokande experiment. 

However, the most natural dark matter candidate in the IDM has a mass moM ^ 60 GeV 
(middle mass range) [9l [11]. Furthermore, being a scalar, it has dominantly spin-independent 
interactions with nuclei. These features open the interesting possibility to enhance the capture 
of the IDM candidate by the Earth, through resonant scattering on heavy nuclei, in particular 
iron. This point has been already emphasized in fT3], an article that appeared while our work 
was being completed. Our results on neutrinos and muons fluxes, obtained using different tools 
but making similar approximations regarding capture in the Earth, are in agreement. However, 
unlike in [13j, we also take into account the constraints from direct detection experiments that 
lead to more pessimistic conclusions. In particular, we compare the flux of neutrinos with the 
expected sensitivity of IceCube, a km^ neutrino telescope under construction at the South Pole, 
and show that this detector may at most give limits that are complementary, but not better than 
those given by existing and forthcoming dark matter direct detection experiments. In an attempt 
to improve this situation, we then consider a scenario based on a natural extension of the IDM 
with three generations of right-handed Majorana neutrinos. We show that, provided we somewhat 
depart from the scheme proposed in [8| for radiative SM neutrinos masses, some scalar candidates 
may have large annihilation rates into mono-energetic SM neutrinos, hence strongly boosting the 
possibility of observing indirect detection of dark matter captured by the Earth. 

Finally, the IDM dark matter may be heavy, moM ^ 500 GeV (high mass range) [TT]. In this 
case, it has a very small cross section on nuclei and is therefore not likely to be captured by the 
Sun or by the Earth. On the other hand, it may accumulate at the centre of the Galaxy and its 
annihilation may produce high energy neutrinos that could be observed on Earth. However, we 
show that the flux is small for mundane astrophysical assumptions regarding the abundance of dark 
matter at the Galactic centre. Here, we set an upper bound on the neutrino flux that might be 
produced by assuming that the annihilation of the IDM candidate does not give gamma rays in 
excess of those observed by the EGRET satellite. 

^This gives us a complete family of WIMPs. One should add the gravitino, a spin 3/2 state, and a spin 2 candidate 
might also exist in theories of modified gravity or as a Kaluza-Klein partner of the graviton in theories with extra 
dimensions. These higher spin DM candidates have essentially gravitational interactions and are not WIMPs. 

^This is a further step toward filling the gap in the table I given in ref. [17j . 
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This paper is organized as follow. The Inert Doublet Model is introduced in section [2j In 
section |3j we fix our conventions on the dark matter abundance in the Galaxy, the Sun and the 
Earth. The section [4] is devoted to a review of capture and annihilation of WIMP dark matter in the 
Sun and in the Earth and their possible signatures in neutrino detectors or telescopes. In section [5} 
we envision the possibility of observing neutrinos produced at the Galactic centre. The predictions 
of the IDM are discussed in section [6) followed by our conclusions. 

2 The Inert Doublet Model 

The Inert Doublet Model is a two Higgs doublet model, Hi and with an unbroken Z2 symmetry 
under which 

Hi Hi and H2 -H2 
and all the other Standard Model particles are even [8l|9]. The potential of the IDM can be written 



as 



V = M?l^lP + M2l^2p + Ai|i/i|^ + A2|i/2|^ + A3|i/lP|i/2p (1) 

+X^\hIH2\^ + ^ [{HlH2f + h.c. 

In this model. Hi contains the standard Brout-Englert-Higgs particle h (the Higgs for short) and 
the discrete symmetry, which prevents FCNC, gives a dark matter candidate in the form of one 
of the neutral components of the extra doublet H2 = {H^, l/y/2{Ho + iAq))^ . The couplings and 
masses of the scalar particles are related by 

ml -2iil = 2\iv^ with v 246 GeV 

^Ifo = + (A3 + A4 + A5)i;V2 

= + (A3 + A4 - A5)i;V2 . (2) 

Depending on quartic couplings, either i^o or can be the lightest Z2-odd particle. We choose 
Ho and, following [9J, we define A^, = (A3 + A4 + A5)/2, the coupling between h and a pair of H^. 
In our investigation of the model we choose 112^ A2 and the masses of scalar particles, including the 
Higgs, as input parameters. 

Experimental constraints on the IDM model are discussed in [9J and further in [21j and [22] . 
The latter also discusses the prospect for discovery of the and H^ at the LHC. In [211, the LEP 
I and II results on the neutralino are used to put constraints on the mass range of i^o and Aq, see 
in particular their figure 8. In the present work, we consider that H^ is the dark matter candidate. 
The masses (and that of H^) we will use are always consistent with LEP data. 

Assuming that i^o were in thermal equilibrium in the early universe, there are essentially 
three distinct Hq mass ranges that are consistent with the WMAP abundance for dark matter. 
In the sequel, we refer to them as the low mass (3 GeV < tuhq ^ 8 GeV), the middle mass 
(40 GeV < rriH^ < 80 GeV), and the high mass (500 GeV < ttiho < 15 TeV) ranges. 

Candidates in the low mass range annihilate only through the Higgs channel. These solutions 
may be compatible with the DAMA result [16j, an hypothesis that may be tested in the future 
by other direct detection experiments (CDMS, XENON and CoGeNT). Candidates in the middle 
mass range annihilate in the Z (when coannihilation with Aq is kinematically allowed) or the Higgs 
channel. Some solutions give a very significant gamma ray line from annihilation at the GC which 
might be observable with the GLAST/Fermi satellite [12j. Above 80 GeV, annihilation into 
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pairs is allowed, the cross section is large and the abundance falls well below the WMAP abundance 
At higher masses, the annihilation cross section tends to decrease and the relic abundance 
increases. This regime is analogous to that of Minimal Dark Matter [23j but, in the IDM case, there 
are more parameters to play with and there is a whole mass range around 1 TeV of candidates with 
an abundance consistent with WMAP [llj. 

3 Dark matter distribution 

The distribution of dark matter in the Galaxy is uncertain, in particular in the central region. 
Rotation curve observations suggest a rather cored profile [241 ES], with a flat behaviour at the 
centre, whereas numerical simulations predict more cuspy profiles in the innermost region of the 
Galactic centre (see Kravtov et al. [25j, Navarro- Frenk- White (NFW) \23\ and Moore et al. [27j as 
well as the recent Via Lactea and AQUARIUS simulations [281 123 )• Other studies even indicate 
the existence of a dark disk [301 ElJ • In the present study, we will work with a fixed astrophysical 
framework and focus on the NFW profile which is parameterized as follows: 



where tq is the distance to the centre of the Galaxy, po = 0.3 GeV/cm^ is the dark matter density 
in the solar neighbourhood and (rg [kpc], ag [kpc], a, /3, 7) = (8.5, 20, 1, 3, 1). In this framework, 
the dark matter halo is spherical and at rest in the Galactic coordinate system while the Sun moves 
on an orbit with mean velocity vq ^ 220 km/s. The velocity of dark matter is supposed to have a 
Maxwell-Boltzmann distribution with mean velocity v = 270 km/s. 

4 Indirect detection from the Sun and the Earth 
4.1 Capture rate 

Dark matter of the Galactic halo might be captured in the Sun or the Earth due to the energy loss 
caused by elastic scattering of dark matter particles on nuclei. The size, density and composition of 
the Sun and the Earth are quite different and so is their ability to capture DM. WIMPs may gener- 
ically couple with a nucleus through axial- vector (spin-dependent, SD) or scalar (spin-independent, 
SI) interactions [32j. In the latter case, elastic scattering with a nucleus is coherent and the cross 
section is proportional to the square of its atomic number A^, while spin-dependent interactions 
depend only on the total spin of the nucleus a oc J (J + 1). The Sun is big and mainly composed 
of hydrogen atoms (^ 77% of the mass of the Sun). Hydrogen is the lightest nucleus and there- 
fore capture is likely to be more important for WIMPs with spin-dependent interactions. On the 
contrary, the Earth is comparatively small and mainly composed of heavy nuclei. Capture of dark 
matter particles may therefore be achieved by scalar spin-independent interactions. 

Our treatment of capture in the Sun and the Earth will be elementary. Both for the sake of 
exploration and because, for the Earth, we will focus on a dark matter candidate in a mass range 
close to the iron mass, we follow the simple approximation of [33] and presented in the classical 
review This approach is expected to work for capture in the Sun. However, the situation for the 
Earth is more complex. The crux of the problem is that the Earth escape velocity is relatively small 
compared to the halo mean velocity, and in general capture in the Earth is sensitive to the (not well 
understood) low velocity part of the dark matter velocity distribution. The approximation of |I33] 
and [3J consists in treating capture of dark matter from the halo as if the Earth was in free space. 




(3) 
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However, the Earth is actually deep in the gravitational field of the Sun. Moreover, either through 
direct collisions or gravitational interactions, dark matter from the halo might be captured on solar 
bound orbits, possibly affecting the abundance and velocity distribution of dark matter particles 
near the Earth position. These and related issues have been investigated in details by Gould in [33] 
and further in [35], by Damour and Krauss [36i|, by Edsjo and Lundberg [37] and, more recently, by 
Peter [38l[39]. To make short a long, albeit fascinating stor}j^- which implies taking into account 
the influence of planets - let us summarize here that these works show that the approximation we 
use for capture in the Earth works reasonably well if the dark matter candidate's mass is close to 
the iron resonance (see figure 16 of [37j and figure 10 of [39j), which is precisely the case we focus on 
in the present work. A more precise approach is clearly possible, for instance using the DarkSusy 
numerical package ^40j, but in the present work we have limited ourself to an analytical treatment 
both for capture in the Sun and in the Earth, an approach which, although not strictly speaking 
correct, is very simple to implement. 

In the plain Inert Doublet Model, elastic scattering of Hq with nuclei only takes place through 
Higgs exchange, which is a scalar, spin-independent interaction, figure 10 The low energy cross 
section for spin-independent scattering of Hq on a nucleus Afi of mass mj\f. is 

The factor /, that parametrizes the Higgs to nucleons coupling, is related to the trace anomaly, 
fmj\f = {J\f \ ^qmqqq\J\f) — ghNN^- This factor is not well known due to the uncertainty on the 
contribution of strange quarks to the nucleon's mass. From the results quoted in ref. [16], we take 
/ = 0.30 as central value for a nucleon (we do not make a distinction between the coupling of the 
Higgs to a proton and a neutron) but it may be within a rather wide range 0.14 </ < 0.66. 



Figure 1: Scattering of on a nucleus M through exchange of a Higgs boson h. 



The rate at which dark matter is captured by spin-independent scattering is given by, 
^0/® _ Plocal 270 km/s sr F f 6 <^ 

^ - 0.3 GeVcm-3 tJ,_, 2^ lO'^o cm^ ' 

where Cq = 4.8 * 10^^ s~-^ (cq = 4.8 * 10-^^ s~-^) sets the scale for WIMP capture by the Sun 
(Earth), ref. [3|. 

The two astrophysical parameters, piocal — 0.3 GeV and viocal — 270 km/s, refer to the local 
density and mean velocity of the dark matter halo. The sum runs over all elements present in 
the respective astrophysical object and encodes the information regarding their composition. The 



^An historical summary, may be found in the article by Edsjo and Lundberg [37) . 

^For candidates in the middle and low mass ranges, inelastic scattering with Z boson exchange (i.e. Hq Aq) is 



severely constrained by direct detection limits which impose ttiaq — rriHo ^ 100 keV 
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parameters and 0^ refer to the mass fraction and the distribution of element i. For those, we use 
the values given in [3j. 

The cross section for capture = cr^(^ifo5 ^A^J spin-independent scattering on nuclei at 
low momentum transfer is given by eq. Q. The Fi = Fi{mHo) are nuclear form factors and the 
suppression factor Si = encodes the kinematics of DM- nuclei elastic collisions. Capture 

in the Earth is greatly enhanced by the appearance of resonances between the masses of the nu- 
clei and of the WIMP. For the discussion, it is useful to introduce the dimensionless parameter 
X = tuho/tuj^.. The effect of resonances can be understood from an interplay of two factors: the 
kinematic suppression Si and the form-factor suppression Fi. The kinematic suppression factor S 
is defined as: 

1 + 

^esc) ^ 1156 km/s {{vfg^) = 13.2 km/s) is the Sun (Earth) escape velocity. 
The behaviour of S depends on the escape velocity and therefore differs for Sun and Earth. S 
asymptotically scales like A when A ^ while 5 ^ 1 for A ^ oc. A is strongly peaked when the 
WIMP's mass matches a nucleus' mass {x ^ 1). For the Sun, the higher escape velocity brings A 
above 1 in all but a few extreme points (small tuhq) for all elements except hydrogen. Because of 
this behaviour, 5 :^ 1 in the biggest part of the mass range. However, in the case of the Earth, 
A ^ 1 almost everywhere except at the resonances. Therefore, 5^1 only in a small region around 
the resonance while elsewhere S ^ A <^ 1. 

The nuclear suppression form factor takes into account the fact that nuclei are made of nucleons. 
If both the WIMP's kinetic energy in the halo and its maximum energy transfer to the nucleus are 
comparable to the characteristic coherence energy, the WIMP sees the substructure of the nucleus 
and the cross section is reduced. For the Earth, it can be derived [41] that this effect is only 
relevant at the resonances and that the suppression is negligible (5% effect) for all elements except 
iron. Therefore, following [3], we assume that the form factor suppression is 1 for all elements but 
iron for which it is parameterized as 

= 1 - 0.26^ . (7) 

In the case of capture by the Sun, one has to take into account that the kinetic energy of the 
WIMP when scattering on a nucleus is not anymore the kinetic energy at infinity (connected with 
the local velocity in the halo) but receives a contribution from the Sun's gravitational potential. 
It can then be shown that the form-factor suppression is no longer confined to the resonances and 
additionally that it applies to all elements j2|: 

= + ii-Fr^)e.p\-('f^r] (s) 

L V logmj, / J 

except for hydrogen, where F = 1. Therein F^^-^ ^ mj, and ai are fit parameters tabulated in [3]. 
The form-factor in the Sun has basically no effect on the lightest elements (hydrogen and helium) 
while for heavier nuclei it leads to a suppression which is increasing with the mass of the element. 
Capture on iron experiences the strongest suppression, 0{W~^ — 10~^) for the heaviest dark matter 
candidates in the middle range under consideration. 

To summarize the dependence on the WIMP's mass and the different nuclei, we define for the 
IDM candidate a function J-'q/^ 

= E/f ^"^^f '^^(— ^^/*(-i^o,Ar.) , "f^ f -g^. (9) 
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The behaviour of J^q/^ as a function of rriHQ is displayed in figure [2] both for the Earth and the 
Sun. It shows in particular the great importance of the iron resonance in the case of the Earth. 
With the definition of ^0/0, the capture rate can be written: 

^0/0 _ ^0/0 Plocal 270 km/s 

^ - ^ 0.3 GeV cm-3 ^^'^^ 




40 60 80 100 

rriHo [GeV] 

Figure 2: Typical behaviour of ^0/0 in function of ttiho for Sun (dashed blue line) and Earth (solid 
black line). The resonance phenomenon, which is present for elements in the Earth, is absent for 
the Sun (here for ruh = 120 GeV, / = 0.3 and Al = 0.2). 



4.2 Annihilation rate 

The number N of WIMP dark matter in the Sun or Earth is controlled by the capture and annihi- 
lation rates 

^ . c - .r, (u) 

where C is the capture rate defined in eq. ([s]) and F is the annihilation rate. For a self conjugate 
particle, like Hq, the total annihilation rate from a (effective) volume Veff is given by 

av. (12) 



2Veff 

The cross section for annihilation of Hq is through the Higgs channel only, figure [3} and is given by 



TT (4m|^ - mi 



Ho - •"h) ^Ho 



With eq. (12) the solution of (11) is 



tanh2(^/— Ctoj = -tanh2(t0/r) = - Feq, (14) 
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Figure 3: The only tree level annihilation channel of light Hq {tuhq ^ '^w) into fermions pairs is 
through the Higgs boson. 



where 



i^EQ = tanh^(t0/T) and r = a/ — (15) 



and ^0 = 4.5 Gyr is the age of the solar system. The effective volume of the Sun or of the Earth 
is determined by roughly matching the core temperature with the gravitational potential energy of 
the WIMP at the core radius. Following [3l ED SI we adopt 

K| = 5.8*10^^ cm^f^)'^', (16) 



VS = 1.8.10- cm3(^)"^ (17) 

If the time scale r for capture and annihilation to equilibrate is small compared to the age of the 
solar system, the equilibrium factor and annihilation rate reach a maximum, ^ r ^ F^q = 

1 ^ r = c/2. 

Given F, the differential flux of annihilation products i in the Sun or Earth is given by 

F ^ /dN 



(i). - 4^E-«.(IL as) 

F 

where R is either the distance from the Sun to the Earth (i? = ^0-0 = 1-5 * 10^ km) or the radius 
of the Earth (i? = r0 = 6300 km). BRp is the branching ratio for annihilation into final state 
F and {dN/dE)F,i is the differential spectrum of the annihilation products i through the decay 
channel F. The sum on F runs over all the kinematically allowed final states. 



4.3 Muon flux 

In the following, only the resulting flux of muonic neutrinos will be taken into account, since they 
provide the best directional information. Electron- or tau-neutrinos produce a particle shower close 
or within the detector volume hence blurring the directional information. 

For neutrinos from dark matter annihilation, one may distinguish two extreme cases referred 
to as soft and hard spectra. The soft spectrum is falling rapidly with the neutrino energy. The 
hard spectrum on contrary is flatter with significant contributions at high energies. Specifically, 
annihilations of Hq into bb quarks give a soft spectrum of neutrinos, while a hard spectrum originates 
either from pairs of VF or Z bosons (if tuhq > raw) or from r^r~ (if mj^Q < mw) channels. 



The detection rate per unit detector area of neutrino-induced, through-going muon events from 
the Sun or the Earth is computed using the charged-current cross sections and the muon range. 
Ignoring detector thresholds, as in reference [3J, the muon flux in the detector is found to be: 

- <Pe/e E E BR,{Nz^),,M, (20) 

i F 

where (fQ = 2.54 * 10~^^ km~^ yr~^ {(p^ = ip^ (^-e /^e)^ = ^0 * 5.6 * 10^) for the Sun 



(Earth), C®/® is the capture rate defined in eq. ^5, 10) and Feq is the equilibrium factor defined 
in eq. ([l5]). The neutrino-scattering coefficients and muon-range coefficients bi are given in [3J as 
ajy = 6.8, ajj = 3.1, bjy = 0.51 and bjj = 0.67. BR^ is the branching ratio for annihilation into 
final state F and {Nz'^)F,i{mHo) is the second moment of the spectrum {dN/dE)F,i of neutrinos of 
type i from final state F scaled by the square of the injection energy Ei^ 

{Nz')F,i{Ei^) - ^ J {^)^,{E,,Ein)EldE.. (21) 

As neutrino detectors have certain thresholds on the minimum neutrino energy detectable, a lower, 
detector dependent, energy bound is implicit in the integration range. 

5 Indirect detection from the Galactic centre 

Simulations indicate that there is a higher dark matter density at the centre of the Milky Way. As 
for the Sun and the Earth, this opens the possibility of indirect detection of dark matter through 
its annihilation products. The expected flux from the Galactic centre depends on the assumed dark 
matter density profile p. 

For a self conjugate dark matter candidate, like the Hq^ the flux of annihilation products i at 
the Earth is 

J = ^roA..,^ ^ Am Ad, (22) 

with j{An) = — \ — [ dn [ ds p^{r{s,e)), (23) 

^OPiocal Jline of sight 

where tq = 8.5 kpc is the distance to the Galactic centre, dN/dE is the energy spectrum of 
annihilation products and av is the total annihilation cross section times velocity. The factor 
1/2 results from the number of ways to have a pair out of n particles ^n^/2j^ The quantity 
J encompasses all the astrophysics and is defined as the integral of p^(r) along the line of sight 
averaged over a solid angle Afi. With a NFW profile, the integration gives for instance J = 1382 
for Afi = 10~^, corresponding to the resolution of EGRET, or J = 536 for = 2.5°, which is the 
resolution of ANTARES, a neutrino telescope which is able to look in the direction of the Galactic 
centre. 

6 IDM results 

Our discussion is divided into the three mass ranges for the Hq candidate. In each case, the relevant 
astrophysical context is distinct. 



^This factor of 1/2 is in agreement with the expression given by [17], |43] and 
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6.1 Low mass candidate: indirect detection from the Sun 



The light dark matter candidate niHo < 10 GeV of the IDM has been studied in [16j in view of 
explaining the DAMA annual modulation signal. Such a candidate has a larger number density 
than heavier candidates. It is however too light to give a large signal from the core of the Earth 
(see figure [2]). Because of its large number density and quite large cross section, annihilations at 
the Galactic centre might be important and, indeed, a large flux of photons is to be expected which 
might be studied by the GLAST/Fermi satellite [1I6I HBJ. At these energies, however, the flux of 
neutrinos is well below the background of atmospheric neutrinos. The only possible relevant source 
of neutrinos might thus be the Sun. 

The result of our analysis, for rrih = 120 GeV, is given in figure [4] where we show both the 



expected flux of muonic neutrinos according to eq. (19) and that of muons from eq. (20) as a 
function of the mass of the dark matter candidate, ttiho^ and the bare mass scale, /X2- The neutrino 
spectra have been computed using WiMPSiM [45, 46j^and are cut at a threshold energy — 2 GeV. 
The fluxes are shown together with the WMAP allowed region in black (which we have computed 
using MICrOMEGAs [47]) for 0.094 < VtoMh^ < 0.129, the current limits set by direct detection 
experiments (in white, from left to right, XENON and CDMS respectively) and the DAMA allowed 
region (in light blue) for / = 0.3 (cf. ref. ^16j, the allowed region remaining after exclusions of all 
other direct searches). 

As emphasized in and [18], the flux of neutrinos produced in the Sun by such a light DM 
candidate might be constrained by the Super-Kamiokande experiment. Its sensitivity is currently 
only set down to ttidm — 18 GeV from studies of through-going muons [49j. It could however be 
extended down to 2 GeV including stopping, partially contained or fully contained muons [20j. In the 
present work, the estimation of the Super-Kamiokande sensitivity (drawn in blue) has been obtained 
horizontally extrapolating the one given in [49j to lower DM masses. This limit has been rescaled 
to a muon threshold of 1 GeV (cf. [50j) which corresponds to the neutrino threshold E^^ — 2 GeV 
assuming an energy loss of ^ 50 % in the neutrino to muon conversion. Compared to the projected 
Super-Kamiokande sensitivity presented in [2D J. our extrapolation is more conservative. 

Since m|^Q — + A^'^^, a light i^o and 112 in the range of both DAMA and WMAP require a 
large, negative coupling within ^ [—0.6, —0.1]. This in turn requires a rather large A2 to insure 
the stability of the potential. Concretely one needs [9J 

A. > Ai / A. . ?(!!%^, ,24) 

for Ai^2 > and A^ < 0. In figure |4j the shaded regions correspond to a coupling A2 > 1, for 
rrih = 120 GeV. In the low mass range, both the annihilation and the capture cross section are 
proportional to A^/m^ oc ^ ^2)5 ^f- equations (Q and (13). When increasing the 



Higgs mass, this ratio must be kept constant not to reduce the neutrino/muon flux. The Super- 
Kamiokande sensitivity is therefore shifted to higher values of ji2 that are increased oc m/^. On the 
other hand, the displacement of the exclusion limit A2 > 1 is only oc y/rrih- The sensitivity lying for 
ruh = 120 GeV in the A2 > 1 region will require even larger values of A2 for higher Higgs masses. 

Figure |4] shows that Super-Kamiokande is potentially able to test the light IDM candidate and a 
part of the DAMA allowed region (from about 3 GeV to 4.8 GeV). However, the required couphngs 
tend to be rather large. 

A similar conclusion is found adapting the analysis of ref. |18] which gives the limit set by Super- 
Kamiokande assuming a DM candidate that annihilates with branching ratio = 1 either into r — r. 



^To take into account the effect of neutrino oscillations, we take the flux of muon neutrinos at the Earth to be 
equal to one third of the total flux of neutrinos produced at the core of the Sun. 
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Figure 4: Expected neutrino (left) and muon (right) fluxes from the Sun, for a neutrino energy 
threshold of 2 GeV. Dark blue line in the muon flux plot: estimate of Super-Kamiokande sensitivity, 
conservatively extrapolated to light WIMPs (see text). 

Colour gradient - log^^g [km~^ yr~^] (left) and log^^g [km~^ yr~^ (right); WMAP area (black 
lines); DAM A allowed region (dashed magenta lines); XENON, CDMS exclusion limits (white lines, 
from left to right); The shaded regions correspond to A2 > 1. 
(Parameters: ruh = 120 GeV, A2 = 1 (dotted) and A2 = vr (shaded), f = 0.3). 



c — c or b — b pairs, the strongest constraint being for annihilations into r leptons. This limit for 
T^T~ rescaled by the actual branching ratio in the frame of the IDM is shown in figure [5] (in blue) 
together with the DAMA allowed region (in light blue) and other direct detection experiments (in 
green) . 

For the Hq candidate, the limit from the tau channel indicates (in agreement with the result of 
figure [4]) that the limit set by Super-Kamiokande is competitive with limits set by direct detection 
experiments (modulo the uncertainties regarding the abundance of dark matter in the Sun versus 
that at the Earth, as usual). 

6.2 Middle mass candidate: indirect detection from the Earth 

This refers to a Hq with a mass typically between 40 and 80 GeV [llj. Since the Hq has essentially 
spin-independent interactions, the best place to look for neutrinos produced in annihilations is 
toward the centre of the Earth. However, direct detection experiments put very stringent limits on 
the SI cross section of WIMPs and so on the potential signatures in neutrino telescopes. 

In this section, we discuss the muon flux resulting from Hq annihilations in the Earth. First 
we consider the Inert Doublet Model in its most minimalistic version and show that some of the 
candidates that have an abundance in agreement with WMAP give a muon flux that almost reaches 
the sensitivity of a km size detector like IceCube. We then consider two extensions. First we discuss 
the possible effect of one-loop corrections, in the spirit of ^12j but conclude that they may not 
significantly improve the signatures of the IDM into neutrinos. Finally, we consider the extension of 
the IDM with three heavy Majorana neutrinos introduced in [8J. This extension opens the possibility 
of annihilations of Hq into two SM u or two u with energy nijy^jy = ttiho, thus strongly (albeit to the 
price of some fine tuning) increasing the potential signature in neutrino telescopes. 
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Figure 5: The region allowed by DAM A (dashed magenta lines) together with limits on the spin- 
independent elastic scattering cross section of Hq on a nucleus from direct detection experiments 
(dash-dotted green lines), cf. [51]. Super-Kamiokande (solid blue line) sets another limit from the 
r^r~ annihilation channel of Hq. This limit results from the one presented in [48] rescaled by the 
IDM branching ratio. 



6.2.1 The minimal scenario 



Applying the results of section 4.3, the muon flux is given by (14). For SI scattering, the Earth 
is expected to have a higher capture rate for dark matter masses near the main resonances. In 
the range of Hq masses under consideration, the dominant resonance is that of iron. The expected 
muon flux is presented in figure [6] for four different Higgs masses. In this figure, we neglect the effect 
of detector thresholds and give the total flux of muons. Hence, this is the maximal signal that one 
may get, of course assuming standard hypotheses regarding the abundance of dark matter at the 
centre of the Earth. The mass splittings for Aq and (8 and 50 GeV respectively) are consistent 
with LEP data [21j. 

The figures show the muon flux (colour gradient) together with the 2a WMAP abundance 
region. The shaded regions are excluded by stability considerations and the LEP I limit on the 
Z width. We give the limit (in white) imposed by direct detection experiments {i.e. CDMS and 
XENON [52J). Finally, the sensitivity of the IceCube detector is shown in blue. It depends on the 
neutrino spectrum, for which a distinction between so-called 'hard' and 'soft' spectra is made. The 
former results from annihilation channels like t^t~, the latter from bb pairs. For the Earth, the 
IceCube best case sensitivity for ttidm = 50 GeV is ^ 15 muons per km^ yr for a hard spectrum 
and ^ 40 muons per km^ yr for a soft spectrum. Here we adopt the best case sensitivity presented 
for example in [50] as a function of m/:)^ for 50 GeV < ttiho < 100 GeV and its extrapolation down 
to candidates with a mass of 30 GeV. Since the dominant annihilation channel is into bb, we use 
the sensitivity for a soft spectrum. 

Annihilation rates and thus muon fluxes are higher for candidates with a mass close to 
This however gives a cosmic abundance too low compared to WMAP. This tension resolves as the 
mass of the Higgs is increased. Our results are in agreement with those of [13]. However, we find 
that the candidates with a potential signature in IceCube correspond to regions of parameters that 
are not favoured since they either give a too small abundance {ruh = 120 GeV and 140 GeV), or 
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Figure 6: Muon flux from the Earth for the minimal scenario (see text). No detector threshold is 
applied. 

Colour gradient - log^^g [km~^ yr~'^]', WMAP area (black lines); XENON, CDMS exclusion limits 
(white lines); IceCube sensitivity soft (blue lines); excluded regions (shaded). 
(Parameters: A2 = 0.2, AttiaoHo = 8 GeV, Attih+Ho = ^0 GeV, f = 0.3^. 



are within the exclusions limits set by current direct detection experiments {ruh = 160 GeV and 
200 GeV). For the sake of completeness, we illustrate the variation of the muon flux obtained by 
varying the Higgs-nucleon coupling form factor / in flgure [7| 

For small values of /, a part of the parameter space with muon flux above the sensitivity is in 
agreement with direct detection. This region is however outside the 2a WMAP area. Relaxing the 
WMAP limit on the abundance (acceptance obtained with the 3a WMAP area) could give a muon 
flux above the IceCube sensitivity for some candidates that are not yet excluded by direct detection 
experiments. 

6.2.2 A non-minimal scenario: effect of loop corrections 

To increase the flux of muons one possible strategy is to envision candidates with a hard(er) neutrino 



spectrum, thus exploiting the fact that the muon flux is oc E"^, eq. (20). The best case is to have 
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Figure 7: Muon flux from the Earth for the minimal scenario for different values of the factor /. 
The figure comes from the down-left plot in figure |6j when varying /. No detector threshold is 
applied. 

Colour gradient - log^^g [^^~^ y^~^]i WMAP area (black lines); XENON, CDMS exclusion limits 
(white lines); IceCube sensitivity soft (blue lines); excluded regions (shaded). 

(Parameters: M = 0.2, AmAoifo = 8 GeV, Attih+Ho = 50 GeV, / = 0.14 (left) and f = 0.66 
(right)). 



annihilations into mono-energetic neutrinos, a solution that we will consider in section 6.2.3 First 
we consider the possibility that Hq may annihilate into 7 — z/z^ at one-loop. 

At one loop the Hq may annihilate into a pair of photons (77) or into a photon and a Z boson 
(7Z). The relevance of these contributions in the IDM model has been emphasized in [12]. In 
turn, the Z may decay into uu. There is another contribution, through a loop with VF-bosons, for 
annihilation into z/z/7. This contribution has not yet been computed in the literature. Here, to 
explore whether it would be worth being computed, we make the reasonable assumption that this 
contribution does not overwhelm that from Z decay. The energy at which the Z boson is produced 
is Ez = TTiHo + m^/ArriHQ if ttihq > Tnz/2 and the spectrum of neutrinos is peaked around Ez/2. 

In |12], the possibility of direct annihilation into photons and Z bosons has been investigated in 
view of observing gamma rays from the Galactic centre and four benchmark models (BM), shown 
here in table [TJ are proposed. 

To compute the expected muon flux for those BMs, all we need is the decay spectrum of a Z 
boson in addition to the information given in table [T] For a Z with an energy Ez, the neutrino 
spectrum is obtained using WiMPSlM. The muon flux for each benchmark model is then calculated 
with eq. (20) using the branching ratios of table [TJ 

The BMs have specific parameters combinations which are optimized to give the right WMAP 
abundance, not specially a large flux of Z and by consequence not a large muon flux from the Earth 
(BM I has a large branching ratio into a Z, but a relatively small annihilation cross section). For 
this, we need, ideally, a large branching ratio into 7Z (like BM I and III) and a large annihilation 
cross section (BM IV and, to lesser extent, BM III). Additionally, a large cross-section via the Higgs 
is required (BM IV and III again) so that capture in the Earth through SI interaction is efficient, 
cf. eq. Q and figure [ij Finally, we want a dark matter mass close to the iron resonance. BM III 
realizes the best compromise, but even in this case the flux is below the sensitivity of IceCube. The 
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signature in gamma rays from the GC, as envisioned in [T2], may benefit from the possibihty of 
introducing a boost factor (BF), taking into account the uncertainty regarding the abundance of 
dark matter at the GC. Although, with a bit of luck, a signal in gamma rays could be seen even 
without a BF, with BF= 10^ — 10"^, the signature would be spectacular. Such a BF may be possible 
at the GC but is unlikely in the Solar system vicinity, and a fortiori in the Earth. 
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5 
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3 


7 


9 


10 


Muon rate 


7 — Z only 


4.12e-06 


1.82e-06 


3.36e-02 


4.75e-03 


[ km~^ yr""*" ] 


7 — Z + tree level 


4.81e-06 


2.63e-05 


3.16e-01 


8.58e-01 



Table 1: Parameters and branching ratios for the four IDM benchmark models presented in [T2] . 
The maximum muon fluxes assuming no detector thresholds are far below the sensitivity of IceCube. 



6.2.3 Another non-minimal scenario: introducing Majorana neutrinos 

A natural extension of the IDM model consists in the addition of right-handed Majorana neutrinos, 
odd under Z2 0. This introduces the possibility to give a mass to the SM neutrinos through 
radiative corrections. The IDM Lagrangian then contains additional terms, 

C D hij {jJi Ho - k H^)Nj + ^MjNjNj + h.c. (25) 

where the index i stands for the SM neutrino generation (i = e, fi, r) and j for the three extra 
Majorana neutrinos {j = 1,2,3). Further discussion about the generation of neutrino masses can 
be found in [8J. In this extension of the IDM, the dark matter candidate may be either a boson, 
the lighter of the two neutral components of the H2 doublet (Ho or Aq), or a fermion, if one of the 
Majorana neutrinos is the lightest odd particle. As in the previous sections, we take the Ho as dark 
matter candidate and focus on the fact that the direct coupling between Ho and SM neutrinos allows 
annihilations of Ho into pairs of SM neutrinos or antineutrinos, with Ej^^j^ = m^fQ, see figure [8j In 
principle, annihilation into vj? is also possible. However, this process is either p-wave suppressed, 
oc 'U^gj, or helicity suppressed, oc m^,. The cross section for annihilation into or the conjugate 
channel, is of the form 

= ir- — o • (26) 

4^ {m\^ + m\Y 
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To simplify our, exploratory level, discussion we focus on muonic neutrinos and suppose that only 
one Majorana neutrino, with mass ttin, enters the cross section. We do not consider the possibility 
of oscillations, which is a good approximation for the Earth. For reference, from here on we take 
rriN = 100 GeV and set the Yukawa coupling to /i = 0.1. These values are chosen so as to boost the 
signal into mono-energetic neutrinos, which requires a rather light Majorana neutrino. At the same 
time the SI cross section has to be large enough to ensure efficient capture in the Earth. Meeting 
all these constraints requires some fine tuning, but is possible for a broad range of parameters. 

Figure 8: Direct annihilation of Hq into neutrinos (or antineutrinos Vi) through exchange of a 
Majorana neutrino Nj leading to a mono-energetic neutrino flux. 



The total muon flux coming only from annihilations into mono-energetic neutrinos is shown in 
figure [9j for different Higgs masses. The plotted IceCube sensitivity curves (blue lines) are the ones 
for a hard neutrino spectrum. The abundance used for the WMAP limit (black lines) has been 
computed with MICrOMEGAs. Note that the shape of the WMAP regions differs from the ones 
of figure |6] as annihilations into neutrinos or antineutrinos are also relevant at freeze-out. 

For all Higgs masses but rrih = 120 GeV, there is a region of the parameter space that is in 
agreement with WMAP and for which the muon flux is above the IceCube sensitivity. As in the 
previous sections, uncertainties in / affect the muon flux. We have set / = 0.3 but spanning over 
the range of / preserves the overlap observed, for example, for rrih = 200 GeV while an increase of 
/ for ruh = 120 GeV brings the WMAP region within the IceCube sensitivity curve. 

Figure 10, which is a scattered plot over parameters of the model, gives the total muon flux from 
the Earth in function of the DM mass. The left panel shows that for the IDM several points above 
the IceCube sensitivity are in agreement with both WMAP and direct detection. Additionally, 
for the sake of comparison, the muon fluxes for the MSSM are shown in the right panel. The 
interesting IDM candidates are peaked around the iron resonance (50 GeV < ttiho < 60 GeV). For 
parameters in agreement with WMAP (dark blue points for the IDM, yellow ones for the MSSM), 
the IDM fluxes are typically about 2 orders of magnitude larger than the MSSM ones and well 
above the IceCube sensitivity. Actually, there are even a few points that are within the sensitivity 
of Super-Kamiokande. 

On the negative side, we have to face the fact that a boost of the signal into mono-energetic 
neutrinos requires a rather light Majorana particle (rriN ^ 100 GeV) and a coupling to SM neutrinos 
of the order /i ^ 0.1. For those values of the parameters, similar to those used in [53j, the one-loop 
contribution to the masses of the SM neutrinos is quite large, rriiy ^ keV range, which clearly 
contradicts the current limits on neutrino masses [54J. One way around is to make use of the fact 
that the absolute sign of fermions is non-observable and to assume that there is another contribution 
to the SM neutrino masses, beyond the IDM (say through non-renormalisable operators) that 
compensates the contribution from the light Majorana neutrino. Albeit this is contrived in our 
case, as it would require a fine tuning of parameters at least at level 10~^, a similar mechanism 
is sometimes invoked, for instance in supersymmetric models with a sneutrino as the dark matter 
candidate [55]. Presumably, we should also make sure that no other dangerous radiative processes 
are important, like ^ e + 7 [53j. This is, however, beyond the scope of the present work. 
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Figure 9: Muon flux resulting from mono-energetic neutrinos from the Earth for the IDM including 
light Majorana neutrinos. 

Colour gradient - log^^g [km~^ yr~'^]', WMAP area (black lines); XENON, CDMS exclusion limits 

(white lines); IceCube sensitivity hard (blue lines); excluded regions (shaded). 

(Parameters: A2 = 0.2, AttiaoHo = 8 GeV, Attih+Ho = ^0 GeV, tun = 100 GeV, h = 0.1, / = 0.3). 



6.3 High mass candidate: indirect detection from the Galactic centre 

For the dark matter mass range under consideration in this section (500 GeV < rriHo ^ 1500 GeV)J^ 
the most promising signal of annihilation products comes from the Galactic centre. Earth and Sun 
can not produce large muon fluxes since the cross section of dark matter on nuclei, which is inversely 
proportional to the square of the mass of the candidate, is too small. 

In the framework of the basic IDM and for this mass range, annihilation into pairs of W or 
Z gauge bosons, or into Higgs boson pairs, are dominant at tree level. Annihilation into fermion- 
antifermion pairs through a Higgs boson, especially into tt, is also possible. However, for the 
parameter range considered here, these channels turn out to be negligible. 

The Higgs boson and the gauge bosons give distinct neutrino spectra. Those originating from 

^Higher masses are in principle possible, but the range given here is conservative, in the sense that we do not 
considered large quartic couplings ^11^ . 
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Figure 10: Scattered Plots of the muon fluxes of IDM and MSSM from the Earth. IDM fluxes result 
from mono-energetic neutrinos from the addition of light Majorana neutrinos. 
left: The blue points (upward triangles) correspond to combinations of parameters that are in agree- 
ment with 2(7 WMAP and direct detection limits, the light blue ones (squares) are in agreement 
with 2a WMAP but the cross section is above the direct detection limit, the orange ones (down- 
ward triangles) are outside the 2a WMAP region but the cross section is in agreement with direct 
detection, the magenta ones (crosses) violate both constraints. The grey points (dots) are excluded 
by other constraints (e.g. vacuum stability). The red line gives the IceCube sensitivity for a hard 
neutrino spectrum. 

right: All plotted points are within the 3a WMAP region. The blue ones correspond to the IDM 
with Majorana extension (dark blue below, light blue above direct detection limit). The green and 
yellow points are predicted by the MSSM (yellow below, green above direct detection limit). The 
hues are different detector thresholds: AMANDA-BIO (1997-1999) hard (sohd) and soft (dashed); 
AMANDA-II (2001-2003) hard (solid) and soft (dashed); IceCube best case hard (sohd) and soft 
(dashed); BAIKAL (1998-2001); BAKSAN (1978-1995); MACRO (1989-1998) hard; SUPER-K 
(1996-2001) soft, see ref. [50] and references therein. 

(Parameters: rrih = 200 GeV, A2 = 0.2, AttiaoHo = 8 GeV, Attih+Ho = 50 GeV, / = 0.3, 
m^v^lOO GeV, h = 0.1). 



W and Z are considered as hard neutrinos: they have a flat spectrum and a significant fraction 
has high energies. Those produced through the decay of the Higgs are comparatively softer. The 
amount of hard neutrinos among the total neutrino flux is primarily determined by the branching 
ratio into VI/+VI/" and ZZ. 

In ref. it has been shown that, for the mass range 500 GeV < tuhq % 1500 GeV, the 

abundance is in agreement with WMAP for tuhq ^ M2- Furthermore, all but a thin band around the 
diagonal is excluded by vacuum stability and perturbativity. In the following, we consider 112 in the 
range mj^Q — 20 GeV ^ /i2 ^ + GeV and parameterize the mHo-/^2 plane by {tuhq] — 
The diagonal in the (rriHo'-, 1^2) representation becomes the horizontal fi2 — tuhq = 0. For the Higgs 
mass and the mass splittings, we follow the choice made in the same reference, ruh = 120 GeV, 
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^^AqHq — 5 GeV and Am^+^Q = 10 GeV. The coupling A2 = 0.2 chosen as in the previous 
sections completes the parameters set. 

The abundance, computed with MICrOMEGAs, is shown in figure [TT| The WMAP region is 
drawn in black, solid (dash-dottet) lines corresponding to the 2a (resp. 3a) range. The exclusion 
limits (shaded region), as well as the WMAP area, lie, as expected, close to the horizonal line, 
/i2 — rriHo = 0. The limit from vacuum stability corresponding to A2 = 0.2 excludes the region 

> ^ifo- exclusion limit below fi2 = '^/fo I'esults from the less strict perturbativity constraint, 
1 < Xi < 47r. 
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Figure 11: Dark matter abundance obtained from MICrOMEGAs. The black lines represent the 
WMAP density range at 2a (solid fines) and 3a (dash-dotted line). The shaded regions are excluded 
by vacuum stability (upper part) and by perturbativity (lower part). 
Colour gradient - log;LO WMAP area (black lines); excluded regions (shaded). 
(Parameters: rrih 120 GeV, X2 0.2, AttiaoHo 5 GeV, Attih+Ho GeV). 



The cross sections for the different annihilation channels, the total annihilation cross section, and 
the branching ratios are also computed with MICrOMEGAs for the given choice of parameters. It 
is found that, over the biggest part of the allowed region, the VF+VF" channel is dominant, followed 
by Z Z. The Higgs channel is only relevant in a region that is excluded by the constraints mentioned 
above, while ti is negligible everywhere. This implies that the neutrino fiux is mainly hard for the 
mass range considered in this section. 

The total fiux of neutrinos from the Galactic centre is computed with eq. ([22]), for the chosen 
combinations of parameters. For this purpose, the integrated spectrum of neutrinos {i.e. the total 
number of neutrinos produced from all annihilation channels) is obtained from MICrOMEGAs. 
This integration can be carried out for different energy thresholds. More precisely, we do not only 
compute the total number of muonic, but also that of electronic and tauonic neutrinos. Given the 
distance from the Galactic centre, after neutrino oscillations, the three generations arrive at the 
Earth with a ratio : z/^ : z/e = 1 : 1 : 1, and we obtain the fiux of muon- neutrinos by averaging 
over all generations. In addition, we use the total annihilation cross section and J taken from 
section [5] for a NFW profile and the opening angle of the ANTARES detector, a neutrino telescope 
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that may look at the Galactic centre. 

The resulting neutrino flux is shown in figure [T2| assuming a detector threshold of 30 GeV, 
characteristic of that expected for ANTARES. In the plot on the left side, the colour gradient is 
representing log^^g [ km~^ yi""^ ] and the black lines correspond to the 2cr WMAP limit. The 
scattered plot on the right side gives the neutrino flux only for models that are not excluded by 
vacuum stability or perturbativity. The colour code corresponds to whether the abundance is in 
agreement with WMAP (dark blue) or not (light blue). The two red lines represent the sensitivity 
of the ANTARES detector found in ref. [56] for the two extremals of a purely hard (solid line) and 
a purely soft (dashed line) neutrino spectrum. 




mj/JGeV] m//o[GeV] 

Figure 12: Expected neutrino flux in the ANTARES detector from the Galactic centre, for a NFW 
profile, and a threshold energy of 30 GeV. 

left: The colour gradient represents log^^g^^^f km~^ J^~'^]- The WMAP region is drawn as black 
lines and the excluded areas are shaded. 

right: Scattered plot of the neutrino flux together with the sensitivity of ANTARES for a soft 
(dashed) and a hard (solid) spectrum. The dark blue points lie within the WMAP area while the 
light blue ones are outside. 

(Parameters: rrih 120 GeV, X2 0.2, AttiaoHo 5 GeV, Atuh+Ho GeV). 

As in the non-excluded area, Hq annihilates dominantly into W^W~ and ZZ gauge bosons, 
with a branching ratio ^ 80%, the neutrino flux is mainly hard and has to be compared to the 
solid sensitivity curve. However, the expected neutrino flux lies between two and three orders of 
magnitude below the expected sensitivity of ANTARES. 

As is usual in this discussions of indirect detection from the Galactic centre, we make use of the 
fact that the dark matter density profile in the innermost region of the Galaxy is poorly known and 
might be steeper than the assumed NFW profile and enhanced by some boost factor. For instance, 
in [12], the predicted gamma flux is boosted by a factor ^ 10^ reflecting those uncertainties. The 
boost may also be partly of particle physics origin. The annihilation of heavy dark matter candidates 
that interact through relatively lighter particles (for instance, the W, the Z and the Higgs in the 
present case) may have a larger cross section for low relative velocities than that computed at tree 
level, a phenomenon called "Sommerfeld enhancement" [53 [58]. From the similarity of the IDM 
with Minimal Dark Matter models, the expected Sommerfeld enhancement in Hq annihilations at 
the Galactic centre is (9(10^) at most but might be non- negligible. Regardless of the possibility 
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of boosting the flux, we must ensure that the signatures are not in contradiction with existing 
observations, say in gamma rays, see for instance ref. [59]. The observations of the EGRET satelhte 
of the Galactic centre, which we take from the left-middle plot of figure 5c in ref. [60j, set an 
upper limit on the gamma rays that may be produced in annihilations of the Hq. Similarly to 
the computation of the neutrino spectrum, we have used MICrOMEGAs to obtain the differential 
photon spectrum, together with J for a NFW profile and an opening angle Afi = 10~^ corresponding 
to the resolution of EGRET (cf. section [5]). For each point in the space of models, we may compute 
the gamma ray flux and compare it to EGRET data to derive a maximal boost factor (for the 
masses considered, the most constraining observations are at an energy 10 GeV). Typically 

the maximal boost factors are ~ 10^ — 10^, which are actually quite moderate values, given the 
uncertainties on the Galactic abundance and the possibility of Sommerfeld enhancement. The 
maximum allowed neutrino flux is obtained according to 



/^EGRET f 

V-'^norm J 



/norm _ /NFW "^7 ^ 



(27) 

for each model. The result is shown in figure [13) In the left plot, the WMAP area is drawn in black 
and the excluded areas are shaded. The scattered plot on the right contains only points that are not 
excluded, together with the soft (given for information) and hard ANTARES sensitivities in red. 
As discussed above, the neutrino flux should be compared to the hard sensitivity (solid line), due to 
the large branching ratio into the VF+VF~ and ZZ channels. The dark blue colour corresponds to 
models that agree within 2cr with WMAP, the light blue ones have too high or too low abundance. 

For i^o candidates with mass > 900 GeV, the maximum muon flux is above the ANTARES sen- 
sitivity and several models have an abundance that is in agreement with WMAP. Similar predictions 
for the neutralino may be found in ref. [59j . 



7 Conclusion 

In this paper, we have investigated the indirect detection signatures in neutrinos, for various 
scalar dark matter candidates within the framework of the Inert Doublet Model. Concretely, we 
have considered three distinct mass ranges. The first one corresponds to a rather light WIMP, 
^ifo ^ GeV. This range is not excluded by existing direct detection experiments, and may 
be consistent with the combined results of DAMA/Nal and DAM A/Libra. We have shown that 
the neutrinos that are produced by a light captured in the Sun can be constrained by Super- 
Kamiokande measurements. This result corroborates the (model independent) statements made in 
refs. [18^ 48j. We have also studied the predictions for more mundane WIMP candidates, corre- 
sponding to tuhq ^ 60 GeV. Since the i7o has only a SI cross section on nuclei, the best hope to 
detect its presence is in the Earth, in particular using the resonant scattering on iron. Such signa- 
tures have also been considered in [13j and our results for the fluxes are in agreement. However, 
direct detection strongly limits the potential of signatures from the Earth (see for instance ref. [61j), 
and the Hq is no exception. In particular, we have found no solutions that give an observable sig- 
nal from the Earth and are both in agreement with the WMAP abundance and direct detection 
exclusion limits. In order to boost the signal in neutrinos from the Earth, we have considered an 
extension with (light) Majorana neutrinos. The model has to be substantially fine-tuned to have 
WMAP abundance, capture in the Earth and strong annihilations into mono-energetic neutrinos 
(as well as reasonable SM neutrino masses), but some solutions may be possible. These are depicted 



in figure 10, Finally, we considered heavy WIMPs, tuhq ^ TeV. Those candidates annihilate at the 
centre of the Galaxy and produce high energy neutrinos that may be observable, provided a boost 
factor - 10^ 
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Figure 13: Upper limit on the expected neutrino flux from the IDM in the ANTARES detector. The 



neutrino flux of figure [12] is enhanced by a boostfactor factor in order to circumvent the astrophysical 
uncertainties on the dark matter profile at the Galactic centre and to obtain an upper limit on the 
expectable neutrino flux while leaving the corresponding 7- flux below the EGRET limit. 
left: The colour gradient represents log^^g km~^ J^~'^]- The WMAP region is drawn as black 

lines and the excluded areas are shaded. 

right: Scattered plot of the neutrino flux together with the sensitivity of ANTARES for a soft 
(dashed) and a hard (solid) spectrum. The dark blue points lie within the WMAP area while the 
light blue ones are outside. 

(Parameters: rrih 120 GeV, X2 0.2, AttiaoHo 5 GeV, Attih+Ho GeV). 



Note added in proof 

Recent simulations of the dark matter halo in presence of a disk of baryons indicate that there might 
be a rotating disk of dark matter in our Galaxy [30]. Because of the small relative velocity between 
this so-called dark disk and the solar system, it has been shown in [62] that capture in the Earth, 
and to some extent in the Sun, might be dramatically enhanced. In that paper, the impact of the 
dark disk on the neutrino and muon fluxes from the Earth is illustrated for CMSSM neutralino 
candidates, but clearly similar enhancements should be expected for IDM dark matter candidates 



{i.e. in particular for a candidate close to the iron resonance, as considered in section 6.2). We keep 
for further work the investigation of this exciting possibility. 
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